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ABSTRACT 

We have re-analyzed archival RXTE data of the X-ray binary Cygnus X-3 with 
a view to investigate the timing properties of the source. As compared to previous 
studies, we use an extensive sample of observations that include all the radio/X-ray 
spectral states that have been categorized in the source recently. In this study we 
identify two additional instances of Quasi-Periodic Oscillations that have centroid 
frequencies in the mHz regime. These events are all associated to a certain extent 
with major radio flaring, that in turn is associated with relativistic jet ejection events. 
We review briefly scenarios whereby the Quasi-Periodic Oscillations may arise. 
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1 INTRODUCTION 

Cygnus X-3 is a well-known X-ray binary (XRB) located 
in the plane of the Galaxy at a distance of approximately 



9 kpc (Dickey fc Lockman 1990 Predehl et al. 2000). Its 



discovery dates back to 1966 (Giacconi et al. 19671 but the 



intrinsic nature of the system remains a mystery despite ex- 
tensive multiwavelength observations throughout the years. 
A truly puzzling feature of Cygnus X-3 is that the 4.8-hour 
periodicity in the X-ray lightcurve ( jParsignault et al.|1972| ) 
attributed to orbital modulation is characteristic of typical 
low-mass XRBs, but infrared observations suggest the mass- 
donating companion to be a Wolf-Rayet (WR) star (van 
Kerkwijk et al.|l992 l, characteristic of high-mass XRBs. 

Unlike most other XRBs, Cygnus X-3 is relatively 
bright in the radio virtually all of the time (> 100 mjy), with 
periods of quenched or quiescent radio emission punctuated 
by giant outbursts of the type unseen in classical radio- 
emitting XRBs (Walt man et al.| [l996). During these out- 
bursts there is strong evidence of relativistic jet-like struc- 



tures ( |Mioduszewski et al.|2001| |Miller- Jones et al.|2004 |. 

Surprisingly, given the wealth of available timing data, 
the X-ray timing properties of Cygnus X-3 are not particu- 
larly well studied. Thus, apart from the strong orbital modu- 
lation present in the X-ray lightcurve, the timing properties 
of Cygnus X-3 appear remarkably nondescript. The power 



density spectrum (PDS) of Cygnus X-3 was studied in Will- 
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ingale et al. ( 1985 1 and was found to be well described by a 
power law of index /3 = — 1.8 in the frequency range 10 -5 - 
0.1 Hz. This result has been more or less verified in further 
studies by Choudhury et al. (2004; f3 ~ —1.5) and Axels- 
son et al. (2009; f3 ~ —2 in the canonical, disk-dominated 
soft X-ray state and —1.8 in the canonical, Comptonization- 
dominated hard X-ray state). The PDS has no power above 
0.1 Hz, a result which was discussed in Berger & van der 
Klis (1994; they placed an upper limit of 12% rms above 1 
Hz) and subsequently in Axelsson et al. (20091. This is be- 



lieved to be an effect of scattering in the nearby surround- 



ing medium (Zdziarski et al. 20101. However, van der Klis 



& Jansen ( 1985 I found short intervals (5-40 cycles) of mHz 



quasi-periodic features in the PDS during the soft state and 
possibly also during the hard state in EXOSAT/ME data. A 
mHz feature in the Cygnus X-3 hard X-ray lightcurve was 
also found in a balloon flight study by Rao et al. (19911. 



However, no Quasi-Periodic Oscillations (QPO) were found 



in the RXTE/PCA study in Axelsson et al. (2009) in ei- 



ther the soft or hard spectral states of the source. It is im- 
portant to note that all of these previously studied PDSs 
are indicative of red noise (the "leakage" of power in low 
frequency features to higher frequencies) dominated PDSs. 
The discrepancies between the reported results clearly show 
that these transient features are short-lived and sporadic or 
time-sensitive events. Interestingly, the previously discov- 
ered QPOs appear to occur following major radio flaring 
episodes ( Koljonen et al.|2011 1. 

In this paper we study the PDS of Cygnus X-3 using 
the extensive archive of RXTE/PCA timing data including 
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pointings during and after major radio flaring events. We 
introduce the data and the analysis method in Section 2 
and present our results in Section 3. In Section 4 we briefly 
review possible causes for the QPOs in Cygnus X-3 and 
summarise the paper in Section 5. 



2 OBSERVATIONS AND METHOD FOR 
CURRENT STUDY 

2.1 Observations 

We analyzed, following standard procedures using FTOOLS 
v6.9, the archived RXTE observations from 1997 to 2011 
totaling 172 pointings and more than 0.7 Msec of data, 
covering all spectral states as defined in Koljonen et al. 
(2010; hereafter K10), and orbital phases. The RXTE/PCA 
lightcurves were extracted using the generic binned data 
with a time resolution of 4 ms from channels 0-35 cor- 
responding approximately to the energy range 2-15 keV 
(epoch dependent). Subsequently, the lightcurves were 
binned to a 100 ms resolution. We verified our results by 
comparing the outcome to Standard 1 data that have a time 
resolution of 125 ms. In order to gauge the X-ray spectral 
stat^Jof the source during the selected pointings (if not al- 
ready done so in K10), we extracted lightcurves from chan- 
nels 3-11 and 24-37 corresponding to energy ranges 3-6 keV 
and 10-15 keV using the Standard 2 data products and com- 
pared the resulting hardnesses and intensities of these bands 
to the hardness-intensity diagram presented in K10. 

Furthermore, in order to identify the correct radio/X- 
ray state simultaneous radio observations were inspected 
using archival data from the Green Bank Interferometer 
(GBI), Ryle/ AMI-LA, and RATAN-600. To follow the over- 
all evolution of the radio/X-ray spectral states within the 
context of the selected pointings containing the QPOs, 
the monitoring data from RXTE/ ASM, Swift/BAT and 
CGiJO/BATSE were used in phase-selected format as in 
K10 in addition to the radio lightcurves from the above- 
mentioned radio observatories. 



2.2 Timing analysis 

The timing analysis in this study was undertaken using ISIS 



(Interactive Spectral Interpretation System, Houck & Deni- 
|cola|2002] ) with the SITAR timing analysis module. The usual 
procedure for searching for quasi-periodic phenomena in 
lightcurves is to construct a power density spectrum (PDS) 
and search for sharp peaks in this spectrum. To construct 
a PDS we prepared the 100 ms time resolution lightcurves 
(done separately for each pointing) and calculated the PDS 
for segments of length 8192 bins over the whole lightcurve 
rejecting the segments with data gaps and finally averaging 
the PDS over the whole lightcurve. To determine the sig- 
nificance level of the possible QPO detections we followed 
the Monte-Carlo analysis of |Benlloch et aL\ ( |2001[ ) which re- 
lies on simulating random red-noise dominated lightcurves 



with the algorithm described in Timmer & Konig ( 1995 1 



First, we performed simple fits to the PDS with ISIS (with 
high degree of binning), where we include a simple power- 
law for modelling the low-frequency noise and a constant 
for modelling the Poisson noise. Then 5000 lightcurves were 
simulated, using the model of the PDS ensuring that the gen- 
erated lightcurves had the same length and mean counts as 
the original observations. Based on the distribution of the re- 
sulting PDSs from simulated lightcurves for each frequency 
bin, 95%, 99% and 99.9% confidence intervals were obtained. 
Also, root-mean square values were obtained from the dy- 
namically calculated PDSs from a frequency range 0.004-0.1 
Hz using the same segment length over the lightcurves as 
above. 

The orbital phase of individual segments of lightcurves 



were determined using a cubic ephemeris ( Singh et al. 2002 1 



3 RESULTS 

3.1 Previous detections of QPOs 

QPOs have been detected at least six times according to the 
literature in the first studies dedicated to this topic (Table 
1). EXOSAT/ME observations exhibited QPOs with ampli- 
tudes between 5 and 20% of the 1-10 keV flux and periods 
in the range 0.7-15 mHz (50-1500 s) which persisted for 5- 



40 cycles and occurred in the phase interval 0.0-0.75 (van 



der Klis & Jansen||1985 1. A balloon study (Rao et al. 111991 \ 



showed a QPO in the 20-100 keV lightcurve with a fre- 
quency of 8 mHz (121 s), pulse fraction 40% and duty cycle 
60% during maximum brightness. 

A later study based on RXTE pointing data from 1996 
through 2000 was conducted by Axelsson et al.| (|2009[) and, 



as mentioned above, they found no evidence of QPOs on any 
timescale. Specifically, on the shorter timescales (> 10 -3 Hz) 
they do not detect QPOs at all and the power density spec- 
trum is well-described by a powerlaw of index —2, while on 
longer timescales (< 10 -3 Hz), they find that the variabil- 
ity is dominated by the state transitions. However, as [Scott] 
et al. ( 2003 I point out in their study of timing noise in the 
Crab, random walk processes will yield a powerlaw index of 
—2. This suggests that the flaring in Cygnus X-3 produces 
the red noise that essentially mimics a random walk process 
and thus hides any potential intrinsic variability. 



3.2 Current study 



1 Following K10 from hard to soft: quiescent, transition, flar- 
ing/hard X-ray (FHXR), flaring/intermediate (FIM), flaring/soft 
X-ray (FSXR) and hypersoft state. 



Our study yielded only two PDSs that show clear peaks 
above the 99.9% confidence limit, as described below. How- 
ever, we would like to note that while the PDSs are consis- 
tent with power law noise with j3 ~ —1.5-2.5 with a mean 
value of —2.0, there are pointings that exhibit more struc- 
ture in the PDS around 0.01 Hz with multiple peaks just 
below the 99.9% significance level. 

To identify the state of the system when these PDSs 
arise, we calculated the fraction of data as a function of 
phase whose root-mean square values are above the level of 
a standard deviation + the median rms value of all the data 
(see Fig. [I] black histogram) . In addition, each fraction was 
further divided among different radio/X-ray states (coloured 
histograms in Fig. [11. Quiescent radio/X-ray states are not 
shown in Fig. [I] as they have lower rms values. Fig. [I] shows 
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Table 1. Log of QPOs identified from Cyg X-3 in previous studies. 



Mission" 


ObsID 


Date 




Phase 


TSMF 6 


Centroid freq. 


X-ray state 






(MJD) 






(days) 


(mHz) 




E 


1983 Oct 22 


45629.11- 


-45629.18 


0.41-0.75 


11.5 


2.6 


high 


E 


1983 Nov 5 


45643.01- 


-45643.15 


0.14-0.54 


25.5 


2/15 


high 


E 


1983 Dec 18 


45686.14- 


-45686.22 


0.13-0.37 


68.6 


1.7 


high 


E 


1984 Jan 3 


45702.55- 


-45602.73 


0.18-0.64 


84.9 


0.7/3.4 


high 


E 


1984 May 21 


45841.64- 


-45841.92 


0.00-0.13 




1 


low 


B 


1986 Mar 18 


46507.1 








8 





a E=EXOSAT (1-10 keV); B=balloon flight (20-100 keV). 
b Time Since Major Flare. "— " corresponds to not known. 
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Figure 1. The fraction of data whose root-mean square values are 
above the level of a standard deviation + the median rms of all the 
data in a bin spanning a tenth of a phase (black histogram). The 
coloured histograms show how the fractions of data are divided 
among states. 



that most of the rms variation is due to the rms-flux relation 
(tracing the orbital modulation), but two anomalous peaks 
show up at phases 0.3-0.4 and 0.6-0.7. The first peak is 
due to only the FHXR state and the second peak results 
from the excesses in the FIM/FSXR states. These states 
and related phases are the same states and phases where 
the QPOs discussed in this work are found. 

Fig. [2] and Table [2] show the PDSs and significances 
(based on the Monte-Carlo analysis referred to in Section 2) 
respectively of the QPOs we have identified in this subset of 
RXTE data with Fig. ^highlighting the data from the 2000 
Apr 3 observation. As noted in Table |2j the QPOs occurred 
during the flaring states - particularly in the FSXR/FHXR 
states. In the context of phase, we note that the QPOs oc- 
curred between phase 0.2-0.7. Interestingly, the QPOs dis- 



cussed in van der Klis & Jansen ( 1985 I all appear exclusively 



in the phase interval 0.0-0.75 that corresponds to the rising 
branch and top of the phase-folded X-ray lightcurve ( |Vilhu| 
et al. 12003} . 




0.01 0.10 0.01 0.10 

Frequency [Hz] 

Figure 2. The panels show the RXTE/PCA PDSs from each 
pointing tabulated in Table[2] Panels a and b show the lightcurves 
from the 2000 Apr 3 pointing and panel c for the 2009 Aug 9 
pointing (binned in 5 sec bins, to marks the MJD and po the 
phase for the first bin) with arrows marking the QPOs. The box 
in panel b shows the area that is zoomed in Figure [3] Panel d 
shows the PDS from the 2000 Apr 3 pointing and panel e from the 
2009 Aug 9 pointing with Monte-Carlo significance and Poisson 
noise levels. See text for details on methods. 



3.2.1 2000 Apr 3 

This pointing occurred 2.2 days after the peak of a major 
flare (~13 Jy in the 15 GHz band) when the source was 
in the FSXR state. We note the presence of one 8.5 mHz 
(~ 120 s) QPO in this pointing, occurring in two separate 
episodes each lasting about 12 cycles (data between these 
two separate episodes are devoid of QPOs). In the PDS this 
QPO corresponds to 3% rms and most of its power is con- 
centrated in the three cycles visible at the beginning of the 
lightcurve in Fig. [2j>. A blow-up of this region is plotted in 
Fig.[3]where there is possible indication of a weaker 30 mHz 
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Figure 3. A blow-up of the beginning of the lightcurve in Fig.|2p 
showing strong pulsations of ~ 120 sec and superposed on that 
two weaker sequences of 30 sec and 10 sec (marked with arrows) 
in alternating fashion. 



(~ 30 s) QPO-like period (seen in the PDS as a multiple- 
peaked region but not producing a peak above the 99.9% 
significance level) and 100 mHz (~ 10 sec) QPO-like pe- 
riod (not seen in the PDS, however possibly contributing 
power only in this short segment of lightcurve) in alternat- 
ing rhythm: three times the 30 sec period which form the 
peak of the 120 sec QPO and four times the 10 sec period 
which form the valley of the 120 sec QPO. 



3.2.2 2009 Aug 9 

This pointing occurred 17 days after a major flare event (~ 
1.5 Jy in the 15 GHz band). During this pointing the 15 
GHz radio flux density is ~ 50 mjy and the source is in the 
FHXR state. There is a strong 21 mHz (~ 50 s) QPO lasting 
about 20 cycles corresponding to 2% rms in the PDS. In the 
PDS there are also two peaks with lower significances that 
correspond to frequencies 9 and 31 mHz. 



4 DISCUSSION 

An in-depth analysis is beyond the scope of this paper. 
Here we speculate briefly on the origin of the QPOs found 
in this study. 

• QPOs and the geometry of the system. If we assume 
that the QPOs are X-ray emitting blobs or regions in 
orbit around the compact object, the possible timescales 
are determined by the innermost stable circular orbit 
(ISCO, highest frequency) and the Roche Lobe radius 
or the bow shock that is formed as the compact object 
plows the stellar wind of the companion (lowest frequency). 
Identifying the Keplerian radius of the QPOs (~O.16R,0) 
with the ISCO results in a compact object of thousands 
of solar masses therefore rendering this scenario rather 
unlikely. At the other end of the length scale, we find 
the Roche Lobe radius to be approximately 0.7-2.1 R0 
and the characteristic length scale of the bow shock, 
Rbs = mv win d/4npVc ~ 3.1-5.7 R©, with the compact 



object moving with velocity Vc through the stellar wind 
of density p. The above values have been calculated using 
the mass ratio of the system (Mwr/Mc = 3.8jlJ' 4 given 
in Vilhu et al. 20091, orbital period of 4.8 hours, compact 
-30 Mq and typical parameters for 
= 1600 km s _1 with 



object mass Mc = 1.4 
wind accretion from WR stars (v wi „d 
velocity law following Springmann 199 4] and mass-loss rate 
ra = 1.6 X 1O _5 M yr -1 following Szostek fc Zdziarski 
20081. Therefore, for a reasonable range of compact object 



masses the QPOs seem not to be associated with these radii. 

• QPOs as disk oscillations? In the literature there 
are several methods as to how low-frequency QPOs can be 
invoked in the accretion disk, e.g. low-frequency "dynamo 
cycles" in the azimuthal magnetic field (O'Neill et al. 
120101) or Accretion- Ejection Instability (Varniere & Tagger 



2002[ ) to name a few. However, the geometry of the system 
might present a challenge to disk-based phenomena. If the 
accretion disk is perpendicular to the jet and if the jet is 
inclined ~ 14 degrees to the line of sight ( Mioduszewski et 
aT pKnt , then we should be viewing the disk face-on. The 



problem then arises as to how the actual modulation can 
be seen at all. A further interesting question is how would 
the simultaneous frequencies observed in the 2000 Apr 3 
pointing arise in the disk, especially when they seem to be 
producing a clear pattern? Are they causally linked and/or 
produced by the same underlying mechanism? Interestingly, 
the 2009 Aug 9 pointing indicates that the peaks in the 
PDS form a 3:2:1 frequency relation, the same relation as 
has been reported from Sgr A* ( |Aschenbach et al.|2004 i. 



• QPOs as coronal oscillations? As the QPOs appear 
to be present in energy bands from 2-15 keV and as they 
appear to arise in flaring radio/X-ray states that can be 
modelled solely with a Comptonization component (K10) 
this could imply that the actual component causing the 
Comptonization of soft photons is oscillating. An oscillating 
corona due to a magneto-acoustic wave propagating within 
the corona producing multiple QPOs has been theoretically 
discussed in Cabanac et al. (20101. 



• QPOs as oscillations due to the jet ? As the QPOs 
appear to occur after major radio flares, the latter most 
likely signaling a jet ejection event, the two are probably 
linked. This connection could be two-fold: either the jet is 
shadowing, therefore modulating, an underlying oscillation 
in or emission from the disk or corona, or the QPOs are 
caused by some structure in the jet (emitting X-rays, e.g. 
a shock). The latter scenario is discussed for blazars in 



Rani et al. (20101 and it could occur in a turbulent region 



behind a shock where dominant eddies with turnover times 
corresponding to the QPO periods could explain the short- 
lived, quasi-periodic fluctuation. It is also possible that the 
jet possesses a helical structure, and a relativistic shock 
propagating down the jet encounters regions of enhanced 
magnetic field and/or electron density intensifying the 
emission and causing fluctuations. 

• QPOs as oscillations due to the wind? Similarly to 
the jet case, clumps in the Wolf-Rayet companion's stellar 
wind could cause a shadowing effect on an underlying oscil- 
lation/emission. The pointings that exhibit QPOs are often 
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Table 2. Log of QPOs identified from Cyg X-3 in this study. 



ObsID 


Date 


MJD 


Phase 


TSMF a (days) 


X-ray state 6 


obs. length (s) 


50062 


2000 Apr 3 


51637.53-51637.62 


0.33-0.78 


2.16 


FSXR 


4526.6 


94328 


2009 Aug 9 


55052.82-55052.84 


0.22-0.34 


17.0 


FHXR 


1974.0 




no. PDS in ave. 


Model p 


Model x 2 /d.o.f. 


Ccntroid freq. (mHz) 


Significance 




50062 


3 


-2.0 


41.7/27 


8.5/30 


>99.9%/~99% 




94328 


2 


-1.8 


22.5/27 


9/21/31 


~99%/>99.9%/~99.9% 





a Time Since Major Flare. 

b Classification of K10 where FSXR/FHXR refer to the "high" state in the canonical nomenclature. 



accompanied by sudden drops in intensity (van der Klis & 



Jansen 1985 their 1983 Nov 5 observation). A quite simi- 
lar 15% drop in the lightcurve can be seen in Fig. [2ja. The 
light crossing time in this drop is about 16 sec and implies 
an upper limit on the size of the emission region radius of 
< 2.4 x 10 11 cm. This value is remarkably close to the sepa- 
ration between the components a ~ 3 x lO n (M/3OM ) 1/3 
cm, and thus could be produced by a light-blocking clump. A 
clumpy wind was established in Szostck & Zdziarski (2008). 



5 SUMMARY 

In this paper, we reported on the detection of QPOs in 
Cygnus X-3. In this initial study, we found two occasions 
in the RXTE data where QPOs were present, and these fol- 
lowed major radio flares. We speculated on various scenarios 
that may explain the QPOs, but a more in-depth study is 
necessary to single out a favoured one. 
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